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ABSTRACT 

We present new results uncovered by a re-analysis of a Chandra observation of the gravitationally 
lensed, low-ionization broad absorption line (LoBAL) quasar H 1413+117. Previous analyses of the 
same Chandra observation led to the detection of a strong, redshifted Fe Ka line from the combined 
spectrum of all images. We show that the redshifted Fe Ka line is only significant in the brighter image 
A. The X-ray flux fraction of image A is larger by a factor of 1.55 ± 0.17 than the optical R-band 
flux fraction, indicating that image A is significantly enhanced in the X-ray band. We also find that the 
Fe Ka line and the continuum are enhanced by different factors. A microlensing event could explain both 
the energy-dependent magnification and the significant detection of Fe Ka line emission in the spectrum 
of image A only. In the context of this interpretation we provide constraints on the spatial extent of the 
inferred scattered continuum and reprocessed Fe Ka line emission regions in a LoBAL quasar. 



1. INTRODUCTION 

Recently, the gravitationally lensed (GL) quasar 
H 1413+117 (the Cloverleaf quasar) was observed for ~ 
40 ks with the Chandra X-ray Observatory (Oshima et 
al. 2001, Gallagher et al. 2002). The original goal of 
this observation (PI: M. Hattori) was to detect a cluster of 
galaxies that is thought to contribute to the lensing of this 
system (Kneib et al. 1998; Chae & Turnshek 1999). X-ray 
emission from the cluster lens was not detected. However, 
a remarkable iron emission feature at i?=6. 21+0.16 keV 
(90% confidence) was discovered in the combined spectrum 
of all lensed images. This redshifted emission feature was 
interpreted as Fe Ka fiuorescence from the far side of the 
quasar outflow. 

H 1413+117 is a gravitationally lensed system consisting 
of four images of a distant (2:=2.55) broad absorption line 
quasar (Hazard et al. 1984; Magain et al. 1988). HST 
NICMOS-2 observations indicate that the lens is located 
approximately at the geometrical center of the four im- 
ages. A flrm spectroscopic redshift of the lens has yet to 
be obtained, however, a cluster of galaxies at a redshift 
of z = 1.7 has been reported to contribute to the lensing 
of this system (Kneib et al. 1998). H 1413+117 was the 
first quasar to be detected in the submillimeter continuum 
and in carbon monoxide emission (Barvainis et al. 1995). 
H 1413+117 is also a low-ionization broad absorption line 
quasar (LoBAL) showing low-ionization lines of Al III and 
Fe HI (Monier et al. 1998). LoBALs are BAL quasars 
that show one or more of the low-ionization lines of Mg 
II, Al III, Fe II and Fe III (e.g., Boroson & Meyers 1992). 
Its first X-ray detection was made with the ROSAT PSPC 
(Chartas et al. 2000). X-ray observations of BAL quasars 
indicate that they are heavily absorbed with measured col- 
umn densities ranging between 10^"'^ and 10^^ cm~^. The 
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extreme X-ray weakness of LoBALs, their low values of 
aox, and their high levels of UV continuum polarization 
are taken to imply that the absorbing gas in the direct 
line of sight to the quasar central engine is Compton-thick 
to X-ray radiation from the quasar nucleus and only scat- 
tered or reflected radiation can reach the observer along 
an indirect line of sight (Green et al. 2001; Gallagher et 
al. 2002). 

Here, we present a re-analysis of the Chandra observa- 
tion of H 1413+117 that places constraints on the spatial 
structure of the X-ray iron line and continuum scattering 
regions of this LoBAL quasar. The main difference be- 
tween our and previous analyses is that we performed a 
spatial and spectral analysis of the individual X-ray im- 
ages whereas previous X-ray studies were based on the 
combined spectrum of all images. Throughout this paper 
we adopt a cosmology in which Hq ~ 75 ^175 km s^^ Mpc~^ 
and go — 0.5. 

2. OBSERVATIONS AND DATA ANALYSIS 

H 1413+117 was observed with Chandra for about 38 ks 
on 2000 April 19. Details of this observation were de- 
scribed in Oshima et al. (2000) and Gallagher et al. 
(2002). For preparing the Chandra data for analysis we 
used the CI AO 2.2 and CALDB 2.12 products provided 
by the Chandra X-ray Center (CXC). To improve the spa- 
tial resolution we removed a ± 0"25 randomization applied 
to the event positions in the CXC processing and employed 
a sub-pixel resolution technique developed by Tsunemi et 
al. (2001) and later improved by Mori et al. (2001). The 
deconvolved Chandra ACIS image of H 1413+117 is shown 
in Figure 1. A comparison between the deconvolved X-ray 
image and the optical HST images of H 1413+117 taken 
from the CfA- Arizona Space Telescope Lens Survey (CAS- 
TLES)^ clearly shows the enhancement of X-ray image A 
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compared to that of the optical image A. The X-ray image 
positions are in agreement within errors with those of the 
HST images indicating that the X-ray emission used in our 
spatial-spectral analysis originates from the lensed quasar 
and not from extended emission of the Icnsing galaxy or 
host galaxy. To account for the recently discovered de- 
cay in the quantum efficiency of ACIS (possibly caused 
by molecular contamination of the ACIS filters), we have 
applied a time-dependent correction implemented in the 
XSPEC model ACISABSl . 1.^ The ACIS quantum efficiency 
decay is insignificant at energies above 1 keV and does not 
affect the main results of our analysis. To estimate the 
X-ray flux ratios we modeled the Chandra images of A, B, 
C and D with point spread functions (PSF's) generated 
by the simulation tool MARX (Wise et al 1997). The X-ray 
event locations were binned with a bin-size of 0"0246. The 
simulated PSFs were fitted to the Chandra data by mini- 
mizing the Cash C-statistic formed between the observed 
and simulated images of H 1413+117. The relative posi- 
tions of the images were fixed to the observed NICMOS 
values taken from CASTLES. We find that the X-ray flux 
ratios in the full (0.2-8keV) band are [B/A]fuu = 0.49 ± 
0.08, [C/A]fuu = 0.35± 0.07, [D/A]/„„ = 0.37 ± 0.07 and 
the number of detected X-ray events in images A, B, C, 
and D were 147 ± 13, 72 ± 10, 52 ± 8, and 54 ± 9, re- 
spectively. These flux ratios are consistent with the values 
obtained from the deconvolved X-ray image shown in Fig- 
ure 1. For comparison the HST R-band flux ratios are 
[B/A]r = 0.94 ± 0.01, [C/A]fi = 0.78 ± 0.01, [B/A]r 
= 0.74 ± 0.01. The R-band magnitudes were taken from 
the CASTLES Web site. The flux fraction of image A, 
[A/(A-hB-hC-FD)], is a factor of 1.55 ± 0.17 larger in the 
full X-ray band than in the HST R-band. 

The earlier spectral analyses of H 1413-f-117 used the 
combined spectrum of all images. In our re-analysis of 
the data we search for possible spectral variations between 
images by extracting a spectriim for image A and a com- 
bined spectrum of images B, C, and D. The low S/N of 
these data did not allow us to extract individual spectra 
for images B, C, and D. The spectrum of image A was 
extracted from a circular region of radius l"ccntered 0"5 
SW from the center of image A to minimize contamination 
from the other images. The combined spectrum of images 
B, C, and D was extracted from a circular region centered 
on H 1413-1-117 with a radius of 2"5, excluding the source 
region of image A. We determined corrections to the an- 
cillary response files by simulating the spectra of point 
sources at the locations of the images within the apertures 
used in our analysis. For our simulations we used XSPEC to 
generate the source spectra and the ray-tracing tool MARX 
to model the dependence of photon scattering with energy. 
The background was extracted from a source-free annulus 
centered on H 1413-1-117 with inner and outer radii of 10" 
and 30", respectively. All errors reported below are at 
the 90% confidence level unless noted otherwise. A vari- 
ety of models were fitted to the spectrum of image A and 
the combined spectrum of images B, C, and D, employing 
the software tool XSPEC vll.2 (Arnaud 1996). Spectral 
fits were performed in the 1-8 keV bandpass. Given the 
relatively low-S/N of the spectra we performed fits using 



the C-statistic (Cash 1979) and compared the results with 
those obtained using the statistic. Spectral fits per- 
formed with the C-statistic do not require binning of the 
data and produce reliable estimates of the best-fit param- 
eters and their errors for low-S/N spectra. One of the 
limitations of the C-statistic is that it cannot provide a 
goodness-of-fit measure. However, simple analytic expres- 
sions do exist for estimating the goodness-of-fit when using 
the statistic. We therefore used fits performed with the 
statistic to evaluate the probabilities of cxcxHxling the 
derived values. For the spectral fits performed with 
the x^ statistic we binned the data to have a minimum 
of twelve counts per bin. We did not find any significant 
differences between the estimated best fit parameters and 
their associated errors provided by the C-statistic and the 
X^ statistic. In Table 1 we list the best fit parameters and 
their errors based on fits using the C-statistic. 

The spectra were initially fitted individually with a 
model consisting of a simple power law with Galactic ab- 
sorption due to neutral gas with a column density of TVh 
= 1.8 X 1020 cm-2 (Stark et al. 1992). The model also 
included neutral intrinsic absorption at 2; = 2.55. The 
fit residuals in image A showed an emission feature be- 
tween 1.65-1.85 keV that led to an unacceptable fit, with 
a x^ = 13.56 for 10 degrees of freedom (dof). Only the 
spectrum of image A shows clear evidence of this strong 
emission line. We modeled the residual feature in image 
A by considering a Gaussian component near 1.7 keV in 
the observed frame. This led to a significant improvement 
in fit quality at the 98% confidence level (according to 
the i^-test, F = 5.6) with a x^ = 4.0 for 7 dof. In Fig- 
ure 2a we show the spectrum of image A. Protassov et 
al. (2002), have recently argued that the i^-test cannot 
be applied when the null values of the additional parame- 
ters fall on the boundary of the allowable parameter space. 
They proposed a Monte Carlo approach to determine the 
distribution of the F statistic. We followed this approach 
and constructed the probability density distribution of the 
F-statistic between spectral fits 1 and 2 of Table 1. Our 
Monte Carlo simulations indicate that the probability of 
obtaining an F value of 5.6 or greater is ^1.5 x 10^^ close 
to the result found by using the analytic expression for the 
distribution of the F-statistic. The best-fit rest-frame en- 
ergy of the emission feature is Enne = 6.3±0.1 keV. Adding 
an unresolved Gaussian line to the combined spectrum of 
images B, C, and D leads to a slight improvement in fit 
quality at the 86% confidence level (according to the in- 
fest). The combined observed-frame spectrum of images 
B, C, and D fit with Galactic absorption, neutral absorp- 
tion at the source, and a power-law model is shown in 
Figure 2b. The best-fit values for the photon indices in 
the spectrum of image A and the combined spectrum of 
images B, C, and D, are T = 1.8±0.8 and F = 1.7 ± 0.7, 
respectively (sec fits 1 and 4 of Table 1). Our fits support 
the presence of an intrinsic absorber towards image A with 
a column density of Nu = (2.3±2.2) x 10^^ cm'^ (see fit 
1 of Table 1). We also find a significant intrinsic absorber 
from the fit to the combined spectrum of images B, C, and 
D of Nu = (31^2) X 10^^ cm~2. The observed column in 
images B, C, and D (see fit 4 of Table 1) may represent the 



^ ACISABS is an XSPEC model contributed to the Chandra users software exchange world wide web-site http://asc.harvard.edu/cgi-gen/cont- 
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column density perpendicular to the outflow rather that 
along the outflow as shown in Figure 3 of Oshima at al. 
(2001). They estimated the column along the flow to be 
Nh ^ 10^^ cm~^. More complex models including ionized 
absorption or partial covering were not considered in our 
analysis due to the low S/N of the available spectrum and 
the high 2; of H 1413+117 which results in several absorp- 
tion edges due to Ne, O, N, and C being redshifted outside 
the Chandra bandpass. 

We performed simultaneous fits to the spectrum of im- 
age A and the combined spectrum of images B, C, and D 
to estimate the significance of the difference in the Fe Ka 
lines. The common spectral model for the fit consisted of 
(a) Galactic absorption due to neutral gas with a column 
density of A^h = 1-8 X 10^*^ cm~^, (b) neutral absorption 
at the redshift of the quasar, (c) a simple power law, and 
(d) a Gaussian component near 1.7 keV in the observed 
frame. The photon index F and the energies and widths of 
the Fe Ka lines were linked as common parameters in this 
spectral fit. The best-fit parameters and 90% confidence 
errors are presented in Table 1. The rest-frame equiva- 
lent widths and 90% confidence levels of the lines in the 
spectrum of image A and the combined spectrum of im- 
ages B, C, and D are EWa = l-27to 5^ keV and EWbcd 
= 0.201q'2o keV, respectively. In Figure 3 we show the 
68.3% and 85% confidence contours of the normalizations 
of the Fe Ka lines versus the rest-frame Fe Ka line energy. 
We chose to show the 85% contours instead of the com- 
monly used 90% contours because the 85% contours touch 
exactly in this case. Based on these contours, the probabil- 
ity that the iron lines in the spectrum of image A and the 
combined spectrum of images B, C, and D have strengths 
that are consistent with each other is less than 0.02. For 
this confidence contour plot we fixed the parameters of the 
power-law component of the model to their best-fit values 
listed in fit 5 of Table 1 . We conclude that the Fe Ka line 
in the spectrum of image A is significantly enhanced by a 
factor of about six from that of the combined spectrum of 
images B, C, and D. 

3. DISCUSSION 

Our re-analysis of the X-ray spectra of the individual 
lensed images of H 1413+117 indicates the presence of a 
microlensing event in image A because (a) the flux fraction 
of image A, [A/(A+B+C+D)], is a factor of 1.55 ± 0.17 
larger in the X-ray band than what was detected in the 
HST R-band, (b) a redshifted narrow iron line is present 
with a relatively high significance only in the spectrum 
of image A. To obtain additional insight on the possible 
energy dependence of the magnification of image A, we es- 
timated relative enhancements of image fiuxes in the con- 
tinuum and iron-line emission components based on fit 5 
of Table 1. Specifically, we find that the ratio of the X-ray 
fluxes of image A and the combined images B, C and D 
in the iron line component is [A/(B+C+D)]fo Ka > 1-7 at 
the 85% confidence level (see Figure 2). The ratio of the 
X-ray flux densities at 1 keV (observed frame) in the con- 
tinuum component is [A/(B+C+D)]cont = 0.8 ± 0.2 (90% 
confidence intervals) . These values differ significantly from 
the HST R-band flux ratio of [A/(B+C+D)]r=0.41±0.01 
which is expected to be less affected by a microlensing 
event produced by a star or group of stars in the fore- 



ground lensing galaxy since the characteristic sizes of the 
X-ray continuum emission regions of quasars are expected 
to be considerably less than the sizes of the optical R-band 
emission regions. Specifically, variability studies indicate 
that the sizes of the X-ray and optical continuum emission 
regions are ~ 4x10^^"^^ cm and ~ 4xl0^^~^^ cm, respec- 
tively (e.g., Chartas et al. 2001; Dai et al 2002; Wyithe 
et al. 2000). We note that since the flux fraction of image 
A in the X-ray continuum is greater than the flux fraction 
in the R band, the large equivalent width of the Fe line 
in image A is not caused by a decrease in the X-ray con- 
tinuum but is more likely produced by an amplification of 
the flux of the line. 

A scenario in which intrinsic variability is the cause of 
the detection of a strong iron line in only image A and a 
significantly weaker line in the combined spectrum of im- 
ages B, C and D requires fine tuning such that the flare 
of the iron line needs to occur within a certain window in 
time for it to be detected with a large EW in only one im- 
age and the flare must last for only a few days. Our simple 
two singular isothermal sphere lens model for H 1413+117 
predicts that image B leads images C, A, and D by 5 days, 
9 days and 16 days respectively. A similar lens model by 
Chae & Turnshek (1999) was found to reproduce the ob- 
served relative amplifications fairly accurately. Another 
strong argument against the intrinsic variability interpre- 
tation is that changes of factors of ~ 6 in the equivalent 
width of a fluorescent Fe Ka line in a Seyfert 1 galaxy 
over periods of a few days have not been detected to date 
(e.g., Markowitz, Edelson & Vaughan 2003; Gliozzi, Sam- 
bruna & Eracleous 2003; Lee et al. 2000; Leighly 1999). 
We anticipate that the amplitude of Fe Ka line variability 
to be less in quasars than in Seyfert 1 galaxies due to the 
larger sizes and luminosities of quasars. X-ray variability, 
at least of the continuum emission in Seyfert 1 galaxies 
and quasars, is found to be anticorrelated with X-ray lu- 
minosity (e.g., Nandra et al. 1997; Turner et al. 1999; 
George et al. 2000; Dai et al. 2003). A more likely ex- 
planation of the observed energy dependent magnification 
in the X-ray band and the enhancement of image A in 
the X-ray band compared to that in the optical is a mi- 
crolensing event in image A. Such an event will not lead 
to time-delayed magnifications in the remaining images 
and therefore can explain the non-detection of the Fe Ka 
line in the remaining images. A characteristic size scale 
commonly used in microlensing analysis is the projected 
Einstein ring radius of the star, (e- Emission regions with 
a size significantly larger than i^^ will not be affected by 
microlensing, whereas, emission regions less than C^e will 
be significantly magnified. For H 1413+117 with lens and 
source redshifts of 1.7 and 2.55 respectively, and assum- 
ing an isolated star, the Einstein-ring radius on the source 
plane is Ce ~ 2 x 10i''(Mstar/Afo)i/2/i-0-5 

The observed enhancement of the Fe Ka line implies 
that the size of this iron line region is Rune ~ C-E- The 
enhanced Fe Ka emission observed in image A may be 
produced by microlensing of Fe Ka fluorescence from the 
far side of the outflow. Assuming the Murray et al. (1995) 
disk-wind geometry for BALQSOs we expect that Fe Ka 
fluorescence from the near-neutral accretion disk, thought 
to be emitted (based on modeling of the Fe Ka line pro- 
files in several Seyferts) from annuli of inner radii of ri„ 
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~ 6 and outer radii of Tout ^ 30 Vg (e.g., Nandra et 
al. 1997), where Vg = GMbh/c^ and Mbh is the mass of 
the black hole of H 1413+117, will be heavily absorbed by 
the quasar outflow. The far side of the wind is a more 
plausible source of the Fc Ka line because our line of 
sight to it does not traverse a significant absorbing col- 
umn. We estimate the size of the microlensed fluorescent 
region of the quasar outflow to be < (e- We have assumed 
a minimum launching radius for the UV BAL absorber of 
r,ni„ - 1 X 10i*^(M6;,/10^1/ey/2 cm based on hydrody- 
namical simulations performed by Proga et al. (2000). To 
estimate Mbh we assumed the lensed bolometric luminos- 
ity ofLsoi = 3.9 X 10''^ /i^i erg s'^ (Granato et al. 1996), 
/i is the magnification factor ranging between 10 and 23 
based on lens models of this system (e.g., Venturini & 
Solomon 2003, Chae et al. 1999), LBoi/LEdd = Lsdci 
is the Eddington luminosity and r] is thought to range be- 
tween 0.1 and 1 in quasars. We find, Lsdd = 3.9 x 10^^ 
/x~^77~^ erg s~^, Vg = 4.5 x 10^^ iJ.~^r]~^ cm, Mbh = 3 
X 10^ At"^T?"^ Mq and rmin= 5.5 x 10^^ /u-^/^r?-^/^ cm. 
We note that with the values we have adopted the prod- 
uct /xry is approximately equal to 1 . The relatively smaller 
magnification of the observed X-ray continuum emission 
compared to the line emission suggests that the size of the 
region producing the observed X-ray continuum emission 
is larger than t^E- 

The processes thought to contribute to the contin- 
uum X-ray emission in radio-quiet quasars include inverse 
Compton scattering of photons in disk coronae by UV pho- 
tons originating in the accretion disk, resulting in an in- 
crease in the energy of the upscattcred photons from the 
UV range to the 2-10 keV range. At rest energies above 
~ 10 keV, Compton scattering of photons from the disk 
corona by the disk becomes significant. Variability stud- 
ies of radio-quiet quasars with similar black hole masses 
to H 1413+117 indicate that the size of the emission re- 
gion responsible for these continuum processes described 
above is of order io^^~^^ cm (e.g., Chartas et al. 2001; 
Dai et al 2002.) This, combined with our conclusion that 
the observed X-ray continuum emission must have a size 
» Ce suggests that the X-ray continuum emission in 
H 1413+117 is not viewed directly but most likely is re- 
flected in a scattering region with a characteristic size of 
rscat » Ce- 

Our constraint on the size of the scattering region is 
supported by recent HST observations of H 1413+117; us- 
ing imaging polarimctry with the HST, Chae et al. (2001) 
detected a 0.07 mag dimming and a change of the linear 
polarization component of image D relative to the other 
images during the June 1999 observation which they in- 
terpreted as microlensing of a portion of a large-scale, 
asymmetric continuum scattering region seen in polarized 
light. They also conclude that the microlensing caus- 
tic during their observations was located away from the 
continuum emitting region. They independently conclude 
that the scattering region is IO^^L^q cm > rgcat > 2 x 
lOi*^(Mstar/M0)^/^^^°'^ cm. Wc also note that based on 
the simulations of Witt, Mao, and Schechter (1995) im- 
ages A and D of H 1413+117 have the largest probabil- 
ity of being microlensed due to the relatively large op- 
tical depths for microlensing events along these lines of 
sight. In Figure 4 we present a sketch of a scenario that 



can explain the observed X-ray spectra of H 1413+117. 
For simplicity the caustic plane is oriented perpendicular 
to the observer's line of sight. In this proposed scenario 
the microlensing caustic is centered near the black hole 
and the observed continuum emission in all images orig- 
inates predominatly from a large scale scattering region. 
We expect the continuum X-ray emission region and the 
Fe Ka fluorescence from the near-neutral accretion disk to 
be significantly absorbed along our line of sight. This is 
supported by the fact that H 1413+117 is a LoBAL quasar 
and recent studies indicate that their X-ray spectra may 
be reflection/scattered dominated (e.g.. Green et al., 2001; 
Gallagher 2002). Based on the anticorrclation between aox 
and 2500 A luminosity (Vignali et al. 2003) for radio-quiet 
quasars (RQQs) we estimate an average value of aox ~ - 
1.6 ± 0.2 for a RQQ with / = = 6.1 x 10^" erg s'^ Hz"! 

2500A ^ 

of H 1413+117. In a recent survey of high redshift lensed 
quasars Dai et al. (2003) calculated aox (oiox corrected) = 
-3.3±0.5(-1.8+g;i^) for H 1413+117 based on the X-ray 
flux of the images B, C, and D (image A was not included 
since it is likely microlensed) confirming the X-ray weak- 
ness of this LoBAL. The estimated values of aox are highly 
dependent on the assumed model for the intrinsic absorp- 
tion. The observed X-ray emission in image A also con- 
tains significantly magnified Fe Ka emission originating 
from the far side of the fiow and less magnified emission 
from a portion of the large scale scattering region. 

Another possiblity is that the microlensing caustic dur- 
ing the Chandra observations was not centered on the 
black hole but located further away than the X-ray con- 
tinuum emission region of the inner accretion disk. In this 
scenario the observed continuum emission in all images 
originates from both the continuum-emitting region near 
the inner accretion disk and the large-scale scattering re- 
gion. The X-ray flux observed in image A also contains 
magnified emission from the far side of the outflow and 
magnified emission from a portion of the large scale scat- 
tering region. Since the microlensing caustic in this case 
does not overlap with the central region of the black hole 
complete obscuration of the emission from the inner re- 
gion of the accretion disk is not required to explain the 
observed energy-dependent magnification. To summarize, 
the various emission regions referred to in this analysis in 
order of smallest to largest regions are: (a) A very small- 
scale (<< (e), possibly obscured Fe Ka fiuorescence re- 
gion, (b) A small-scale ($C-e) portion of the reflection 
region from the far side of the outflow, which is relatively 
cold and neutral, slightly redshifted, and is highly ampli- 
fied by microlensing. This region causes the large equiva- 
lent width for the iron line in image A. (c) A small-scale 
( < Ce) X-ray continuum emission region originating near 
the inner-accretion disk, (d) A mid-scale (^ Ce) optical 
continuum R band emission region, and, (e) A large-scale 
{> Ce) X-ray and optical scattering region which is highly 
ionized and reflects the X-ray source without creating a 
fluorescent iron line. 

Optical spectra of H 1413+117 taken with the SILFID 
spectrograph at the Canada- France-Hawaii (CFH) tele- 
scope in 1989 (Angonin et al. 1990) and the HST Space 
Telescope Imaging Spectrograph in 2000 (spectra taken by 
E. M. Monier cited in Chae et al. 2001) indicate that a 
microlensing event along the sight-line of image D was in 
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progress for at least a duration of 11 years. This timescale 
is consistent with the estimate by Witt, Mao, & Schechter 
(1995) of the time for the source to cross an Einstein radius 
for H 1413+117. We may, therefore, expect that the pos- 
sible microlensing event in image A persists for a similar 
timescale. 

The detection of the iron line is possible perhaps be- 
cause of an ongoing microlensing event in image A. Future 
observations of H 1413+117 with Chandra may provide 
a confirmation of the microlensing interpretation if a sig- 
nificant variation in the continuum and/or iron line flux 
ratios in the images is detected. This confirmation would 
validate our microlensing models and our conclusions re- 
garding the relative size of the continuum and iron line 
scatterers. 

We acknowledge financial support from NASA grants 
NAG5-9949, NAS 8-38252 and NAS 8-01128. We thank 
the anonymous referee for helpful comments and sugges- 
tions. 
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CHARTAS ET AL. 



TABLE 1 

RESULTS FROM FITS TO THE SPECTRA OF H 1413+117 



Fit 
1 



Source^ 

H 1413+117^ 
(Image A) 



H 1413+117 



Model 

PL, neutral absorption 

at source, and a Gaussian 
emission line at source. 
The fit is performed to the 
spectrum of image A of H1413. 



PL, neutral absorption 



Parameter 
JVh.a 

E,ine 
feline 

C-statistic / nbins^ 



Value'' 



(Image A) 


at source. 


AfH,A 




The fit is performed to 


C-statistic / nbins 




the spectrum of image A of H1413. 










H 1413+117 


PL, neutral absorption 


^BCD 


(Images B + C + D ) 


at source, and a narrow Gaussian 


^H.BCD 




emission line at source. 






The fit is performed to the 






the combined spectrum 


EWH„e 




of images B, C, and D. 


C-statistic/nbins 














H 1413+117 


PL, neutral absorption 


^BCD 


(Images B + C + D) 


at source. 


-/Vh,bcd 




The fit is performed to 


C-statistic / nbins 




the combined spectrum 


xV"- 




of images B, C, and D. 


PixV'^) 



g+0.8 

(2.3±2.2) X 10^^ cn 
6.3to \ keV 
< 0.14 keV 
l.ltli keV 
4.5/13 
4.0/7 
0.78 

r, 7 + 0.7 

-0.7 

(5±2) X lO^-'' cm~ 
19.2/13 
13.6/10 
0.19 

1.7 f 
(Stl) X 10^^ cm" 
6At°o \ keV 
0.01 keV f 

0.35iS:35 keV 
10.64/12 
7.42/8 
0.49 



(4: 



, 7+0.7 
' -0.7 

X 10--'' cm" 
" 10.64/12 
9.87/9 
0.36 



H 1413+117 
Simultaneous fit 
to images A and 
B+C+D 



PL, neutral absorption 

at source and a narrow Gaussian 

emission line at source. 

The fit is performed 

simultaneously to the 

the combined spectrum 

of image A and 

the combined spectrum 

of images B, C, and D. 



r 

JVh 

EW^ 

EWbcd 
C-statistic / nbins 

xV- 
PixVt^) 



1.6 



+0.8 
0.9 

10^- 



6.3j:° J keV 
0.06i°;?g keV 



1.27 
0.20 



keV 



+ 0.59 
-0.59 
+ 0.36 . y 
-0.20 

22.5/27 
18.0/27 
0.52 



^AU model fits include Galactic absorption towards the source (Stark et 
al. 1992). 

^All absorption-line parameters are calculated in the quasar rest frame. 
^All errors are for 90% confidence with all parameters taken to be of 
interest except absolute normalization. All best fit parameters and as- 
sociated errors are taken from fitfi using the C-statistic. 
"^For results on the combined Hpeclrum of all images of H 1413+117 see 
also Oshima et al. (2000) and Gallagher et al. (2002). 

nbins is the number of pulse height analysis (PHA) bins in the spec- 
trum. 

^ P{x^ /^) is the probability of exceeding for i/ degrees of freedom. 
^ The symbol f indicates that the parameter is frozen in the spectral fit. 
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Fig. 1. — The deconvolved X-ray image of the Cloverleaf in the 1-8 keV bandpass. To improve the spatial resolution we employed a 
sub-pixel resolution technique developed by Tsunemi et al. (2001) and later improved by Mori et al. (2001). The X-ray flux fraction of image 
A, [A/(A+B-f C-l-D)], is larger by a factor of 1.55 ± 0.17 than the optical R-band flux fraction, indicating that A is significantly enhanced in 
the X-ray band. 
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Fig. 2. — (a) The top panel shows the Chandra observed-frame spectrum of image A of H 1413+117 fit with Galactic absorption, 
neutral absorption at the source, a power-law continuum and a Gaussian emission line. In the lower panel the residuals of fit 
1 of Table 1 indicate that this model can account for the spectral features in H 1413-1-117. (b) The top panel shows the Chandra 
observed-frame spectrum of the combined images B, C, and D of H 1413+117 fit with Galactic absorption, neutral absorption 
at the source, and a power- law continuum. In the lower panel we show the residuals of fit 4 of Table 1. The spectra have 
been binned for illustration purposes only. The significance of the Fe Ka line in each spectrum was inferred from fits using the 
C-statistic that does not require binning of the data (see §2. for more details). 
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Fig. 3. — 68.3% and 85% confidence contours of the normalizations of the Fe Ka lines in the spectrum of image A (solid contours) and 
the combined spectrum of images B, C, and D (dotted contours) versus the rest-frame Fe Ka line energy. The parameters of the power-law 
component of the model were fixed to their best- fit values listed in fit 5 of Table 1. 
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Fig. 4. — A proposed representation of a BALQSO disk wind overlayed with a rriicrolensing caustic (dashed hne). This model 
represents the first scenario in tlie discussion section. Tlie caustic plane is oriented perpendicular to the observer's line of sight. 
In our proposed model Fe Ka emission produced from reprocessing off the BAL wind is located within the magnification caustic 
and suffers a significant magnification. Continuum X-ray emission produced from scattering off gas-clouds and the BAL wind is 
extended beyond the diamond-shaped caustic and suffers less magnification. 



